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ABSTRACT 

We study properties of dark matter halos at high redshifts z = 2 — 10 for a vast range of masses 
with the emphasis on dwarf halos with masses 10^ - lO^/i-^M©. We find that the density profiles 
of relaxed dwarf halos are well fitted by the NFW profile and do not have cores. We compute 
the halo mass function and the halo spin parameter distribution and find that the former is very 
well reproduced by the Sheth & Tormen model while the latter is well fitted by a lognormal 
distribution with Aq — 0.042 and a\ = 0.63. We estimate the distribution of concentrations for 
halos in mass range that covers six orders of magnitude from 10^ H^^Mq to 10^"^ /i^^Mq, and 
find that the data are well reproduced by the model of Bullock et al. The extrapolation of our 
results to z = predicts that present-day isolated dwarf halos should have a very large median 
concentration of ~ 35. We measure the subhalo circular velocity functions for halos with masses 
that range from 4.6 X IO^/i-^Mq to 10" /i-^Mq and find that they are similar when normalized 
to the circular velocity of the parent halo. Dwarf halos studied in this paper are many orders 
of magnitude smaller than well-studied cluster- and Milky Way-sized halos. Yet, in all respects 
the dwarfs are just down-scaled versions of the large halos. They are cuspy and, as expected, 
more concentrated. They have the same spin parameter distribution and follow the same mass 
function that was measured for large halos. 



Subject headings: cosmology: theory - 
— methods: numerical 

Introduction 



■ cosmology: dark matter — galaxies: formation — galaxies: halos 



Dwarf galaxies with virial masses 10^ — lO^M© 
are the smallest virialized objects in the universe 
that show evidence of dark matter (Mateo 1998, 
and references therein) . Some of these galaxies are 
so small that they can be detected only in the Lo- 
cal Group. Large amounts of dark matter along 
with the proximity to the Milky Way make dwarf 
galaxies an ideal laboratory for testing the hierar- 
chical cosmologies such as the standard ACDM. 
The rotation curves of dwarf irregular galaxies 



and the velocity dispersions of dwarf spheroidal 
galaxies have been used in the last decade to con- 
strain the density of the dark matter in the central 
parts of these galaxies and, thus, to constrain the 
structure of dark matter (DM) halos that host the 
galaxies. So far the results indicate that DM ha- 
los should have flat cores (Carignan & Freeman 
1988; Carignan & Beauliu 1989; Flores & Pri- 
mack 1994; Moore 1994; Burkert 1995). This 
is inconsistent with the cold dark matter theory, 
which predicts cuspy cores with p oc r~^ (Dubin- 
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ski & Carlberg 1991; Navarro, Frcnk, & White 
1995, 1996, 1997, hereafter NFW) or p oc r'^-^ 
(Moore et al. 1999). The controversy of the flat 
versus cuspy DM cores is still not resolved. Recent 
observations and analysis of dwarf and low surface 
brightness (LSB) galaxies rotation curves continue 
to suggest a dark matter halo density profile with a 
relatively flat core (de Blok, Bosma, & McGaugh 
2003; Weldrake, de Blok, & Walter 2003), but 
sec also S waters ct al. (2003); Rhee, Klypin, & 
Valenzuela (2003). 

A second somewhat related problem faced by 
the CDM cosmogony is the excessive substructure 
predicted by this model in Milky Way-sized ha- 
los as compared with what is observed (Klypin 
et al. 1999b; Moore et al. 1999). Within 
the CDM framework this problem finds a nat- 
ural explanation in the reionization of the uni- 
verse: satellites form only inside subhalos that col- 
lapse early, before the universe was fully ionized 
(Bullock, Kravtsov, & Weinberg 2000; Somerville 
2002; Benson et al. 2002). There may be a differ- 
ent reading of the substructure problem which de- 
pends on how peak circular velocities, Wmax, are as- 
signed to dwarf spheroidal (dSph) galaxies (these 
galaxies comprise most of the Milky Way satel- 
lites). If these galaxies have a larger Wmax than 
what was previously expected, there is no dis- 
agreement: the small number of satellites would be 
explained by the relatively scarcity by which large 
subhalos are formed in the CDM model. In this 
case every large DM satellite hosts a dwarf galaxy. 
The absence of a stellar component in the nu- 
merous small DM satellites could be explained by 
feedback and/or reionization. The problem seems 
to be degenerate: observed dwarf galaxies can be 
hosted by either a relatively low-concentration but 
high-tijnax halo (Hayashi et al. 2003; Stoehr et 
al. 2002) or by a high-concentration and low- 
Vmax halo (Lokas 2002). Thus, the concentra- 
tion of dwarf DM halos is quite important. Accu- 
rate measurements of the concentration for a large 
sample of dwarf DM halos is one of the goals of this 
paper. 

Attempts have been made to simulate dwarf ha- 
los with high resolution within a cosmological con- 
text. Moore et al. (2001) simulated one "Draco"- 
sized halo (Mvir ~ 10^h~^MQ) with several mil- 
lion particles and with high force resolution. Their 
fit of the density profile favor the inner slope —1.3. 



Ricotti (2002) presented several simulations in- 
cluding one with a small (1 /i~^Mpc) box. The 
dwarf halos in this small box are resolved with 
only few tens of thousands of particles at high 
redshifts. Ricotti (2002) finds that at z ~ 10 
the density profiles have inner slopes in the range 
[-0.4,-0.5]. Recently, Cen et al. (2004) repro- 
duced results of Ricotti (2002) using halos with 
very small number (few thousands) of particles. 
Cusps in these simulations are much flatter than 
that quoted by Moore et al. Navarro et al. (2003) 
made high resolution simulations of dwarf halos, 
but those were presented only at z = 0. Dwarf 
halos were found to have steep proflles. Thus, 
the issue of the evolution of the cusps of dwarfs 
is still not settled. Numerical effects are of special 
concern for low resolution simulations. Yet, there 
are other issues, which require attension includ- 
ing treatment of non-equilibrium features in the 
density distribution. 

In this paper, we also study cusps of dwarf ha- 
los. But cusps is only one of the aspects, which we 
are interested in. The statistical properties of this 
population of halos such as the mass function or 
the spin parameter distribution are also studied. 

Dwarf halos with masses Mvir ~ 10''^—10^h~^MQ 
studied in this paper form at very high redshifts. 
If a halo of this low mass remains isolated, it is 
expected to grow only 2-3 times since redshift 
z 3 until z = (van den Bosch 2002) with most 
of the mass ending up in the outer part of the 
halo. Thus, the halo profile should not change 
much since z = 3. The dwarfs may be accreted 
by larger halos and become satellites. If that hap- 
pens, they are significantly stripped and lose most 
of their mass in outer regions. In any case, it is 
important to know what was the structure of the 
dwarfs at high redshifts. This is why we focus on 
the structure of the dwarfs at z ~ 3. 

The paper is organized as follows. In §2 we 

present the cosmological model and describe nu- 
merical simulations used in our analysis. The mass 
function of simulated halos span a mass range of 
seven orders of magnitude. Density profiles of the 
most resolved dwarf halos from the 1 h~^Mpc box 
are presented in §3.1. The concentration versus 
mass diagram along with the analytical prediction 
by Bullock et al. (2001a) is also shown there. In 
§3.2 we compute the subhalo velocity function for 
our most massive halos and show the almost self- 
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Table 1 
Parameters of simulations 



Box 


Mass resolution 


Force resolution 


(78 


Name 


(/i-iMpc) 


(ft-iMo) 


(comoving /i~^pc) 






1 


4.9 X 10^ 


61 


0.75 


A 


25 


1.2 X 10^ 


191 


0.90 


B 


80 


4.0 X 10^ 


610 


0.90 


C 


60 


1.1 X 10^ 


2000 


1.00 


D 



similar nature of the dark matter inside virialized 
systems, in agreement with the results of Moore et 
al. (2001) and Dc Lucia ct al. (2003). The differ- 
ential mass function of halos from the 1 hr^^lpc 
box is calculated in §4.1. A comparison is made 
with the prediction by Shcth & Tormcn (1999). In 
§4.2 we calculate the spin parameter distribution 
of dwarf halos for the 1 /i~^Mpc box and compare 
it with one at z = with the same cosmology. Wc 
discuss some of the results presented in previous 
sections and present our concluding remarks in §5. 

2. Numerical Simulations 

We use a series of simulations of a low-density 
flat ACDM cosmological model with the following 
parameters: fio = 0.3, fix = 0.7, and h = 0.7. 
All simulations are done with the Adaptive Re- 
finement Tree (ART) code (Kravtsov et al. 1997). 
The ART code achieves high spatial rcsohition by 
refining the base uniform grid in all high-density 
regions with an automated refinement algorithm. 

For our code a convergence study shows that 
density deviations less than 10% are expected at 
radii larger than 4 times the formal force reso- 
lution or the radius that contains 200 particles, 
whichever is larger (Klypin et al. 2001). This is 
recently confirmed in a detailed analysis of halo 
profiles by Tasitsiomi ct al. (2003). With the aim 
of studying the structure of dwarf halos, we per- 
form a simulation of an 1 /i~^Mpc box on a side 
with 256'^ particles. One of the motivations for 
choosing this small box is to reproduce the same 
configuration as in Ricotti (2002). The simula- 
tion has been stopped at z = 2.3 when the longest 
fiuctuations in the box are still well in the lin- 
ear regime. At redshift z = 3.3 in the simula- 
tion there are 4 large halos resolved with more 



than 200, 000 particles and with proper force res- 
olution 14/i^^pc. For analysis we also use many 
more smaller halos. 

To complement our study we use three more 
simulations. We use the simulation presented in 
Klypin ct al. (2001), who focused on a con- 
vergence study of the density structure of three 
Milky Way-sized halos at low redshifts. At 2; = 3 
the three largest halos in that simulation have 
a mass of about 10^^ H'^Mq and are resolved 
with ~ 10^ particles. For our analysis we also 
use 12 large halos with more than 5000 particles 
(M > 6 X 10^ h-^Ue). 

We also use a 80 h~^Mpc box simulation, in 
which at redshift z = 3 wc identified 5 large halos 
resolved with more than 10^ particles. At z = 
the halos have become a cluster-sized halo of virial 
mass Mr.,r « 2 X lO^/jr^MQ. The high redshift 
halos are used for our analysis. 

Finally, we use the 60 /i^^Mpc simulation with 
256'^ particles described in Colin et al. (1999); 
Klypin et al. (1999a). The parameters of all the 
simulations are presented in Table 1. In column 1 
we show the comoving size of the computational 
box. The mass resolution is given in column 2, 
while the formal spatial resolution (size of a cell 
in the finest refinement grid in comoving units) is 
shown in column 3. Column 4 gives the value of 
CTg (the rms of mass fluctuations estimated with 
the top-hat window of radius 8/i~^Mpc). Finally, 
column 5 gives the name of the simulation (which 
is also used for objects studied in this simulation). 

In our simulations the halos are identified by 
the Bounded-Density-Maxima (BDM) algorithm 
(Klypin & Holtzman 1997; Klypin et al. 1999a). 
The BDM algorithm first finds positions of lo- 
cal maxima in the density field. These density 
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maxima are found using a top-hat filter with a 
"search radius" which is several times larger than 
the force resolution. Once centers of potential ha- 
los arc found, the algorithm identifies halos around 
them and removes particles which are not bound 
to those halos. This procedure also detects sub- 
halos of larger objects - halos inside halos (for ex- 
ample, satellites of galaxies or galaxies in clusters). 
Particles of a subhalo are bound to both the sub- 
halo and to the larger halo. 

3. Halo structure 

3.1. Density profiles and concentrations 

We start our analysis of halo profiles by study- 
ing all most massive halos (subhalos are not in- 
cluded). The profiles are not averaged over time 
and halos are not selected to be quite or not 
to have large merging events. Thus, some non- 
equilibrium features are expected. Figure 1 shows 
the density profiles of the ten most massive halos 
in the l/i~^Mpc box at z = 3.3 plotted from the 
radii which contain at least 200 particles to their 
virial radii. The profiles are given in normalized 
units: the radius is measured in luiits of the scale 
radius rg, the radius where the logarithmic deriva- 
tive of the NFW profile (see equation 1) is equal 
to —2, whereas the density is measured in units 
of a characteristic density pg = Po/4, the NFW 
density evaluated at r = . These parameters are 
taken from their corresponding NFW fits. The ha- 
los show a variety of profiles. Some profiles have 
significant deviations from the NFW fit (e.g., the 
second and the forth from the top) while; others 
are fitted by NFW reasonably well (e.g., the fifth 
halo). 

The wiggles in the profiles seen in the outskirts 

of most halos can be explained as substructure 
that have not yet reached equilibrium with the rest 
of the halo or as the result of a significant merger. 
For example, halos number 2 and 9 from top are in 
the process of a major merger. Both halos show a 
density enhancement in the peripheral part (more 
pronounced in case of halo 2) . The almost power- 
law profile of halo 9 may be a consequence of the 
way the merger has developed ^. 

^Out of the fifteen biggest halos in the 1 /i~^Mpc box, the 
merger companion of halo 9 has the most peculiar density 
profile. It has the highest value even though it is the 
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Fig. f. — Density profiles of the ten most massive 
halos of the 1 ft,~^Mpc box in normalized units. 
The dotted curves iare the NFW profile. For clar- 
ity we shift down the profiles by —0.5 in logarithm. 
Halos are ordered by mass with the most massive 
halo being at the top of the plot. The density 
profiles are cut at the virial radii. The innermost 
bin contains more than 200 particles. No sign of a 
shallow core is seen in any of the halos. However, 
most halos have wiggles in the peripheral regions, 
which are due to residual substructure and ongo- 
ing merging. The labels An, with n = 1, 4, rep- 
resent the "equilibrium" halos used for subsequent 
analysis. 

No matter how large or small are the deviations, 
the profiles do not show a sign of flat cores in the 
central part of dwarf halos. In other words, steep 
central cusps of halos in our simulations contra- 
dict the results of Ricotti (2002), who finds that 
central slopes of halos in his 1 /i~^Mpc simulation 
are very shallow: p (x r", a ~ —0.2. 

There are some differences in the way how we 
and Ricotti (2002) approximate profiles. Ri- 
cotti does the approximations at much earlier time 
2 10 when halos are much smaller and are less 



second least massive halo on the list. 
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Table 2 

Parameters of the most massive halos at z sa 3 





(/i-ikpc) 


AT 


^max 

(km/s) 


rs 

(proper h~ 


ikpc) 


Po 

(proper MqH'^ pc^'^) 


Halo name 


4.6 X 10^ 


9.3 


933636 


48 




2.0 


5.0 X 10-2 


Al 


1.4 X 10^ 


6.6 


285651 


34 




1.1 


7.8 X 10-2 


A2 


1.1 X 10^ 


6.1 


227275 


28 




1.1 


6.6 X 10-2 


A3 


1.0 X 10^ 


6.0 


208005 


31 




1.0 


7.6 X 10-2 


A4 


2.9 X 10" 


41.7 


238083 


190 




7.1 


5.8 X 10-2 


Bl 


1.2 X 10" 


30.8 


95899 


128 




12.2 


9.2 X 10-3 


B2 


1.1 X 10" 


29.9 


86832 


136 




4.4 


8.2 X 10-2 


B3 


1.2 X 10^3 


140.0 


288061 


647 




26.9 


4.8 X 10-2 


CI 


6.9 X 10^2 


117.3 


175882 


497 




57.2 


6.0 X 10-3 


C2 



resolved. Instead of making two-parameter fits as 
we do, Ricotti uses a very simple approach: he 
uses "location and value of the maximum circular 
velocity" . When we apply the same prescription 
to our high resolution profiles at z = 3, we get 
very noisy results with fits going either above or 
below actual profiles. Yet, we did not get sytem- 
atic effects. In Figure 2 we show the density pro- 
files of the 6 most massive halos at z — 10.5 fitted 
with NFW profiles using two free parameters. The 
NFW fits approximate reasonable well the density 
profiles. Even without any fits it is clear that none 
of the profiles are even close to the shallow profiles 
suggested by Ricotti. 

In the subsequent analysis we mostly focus on 
equilibrium halos, which are fitted by the NFW 
profile. Table 2 presents parameters of the most 
massive halos identified in the simulations A - C 
at redshift z ~ 3. In column 1 we show the mass 
Mvir within the virial radius (column 2). Here 
and below we present radii and densities in proper 
units. The virial radius, Tvir, is defined as the 
radius where the average halo density is S times 
the background density according to the spherical 
top-hat model. Here (5 is a number that depends 
on epoch and cosmological parameters (I^Oi^^a); 
for a flat ACDM model, 5 ~ 180 at z = 3. In 
column 3 the corresponding number of particles 
is given. The maximum circular velocity Umax 

1/2 

= (GM(< r)/r)^^^ , where G is the gravitational 
constant and M(< r) is the mass within the radius 
r, is located in column 4. In column 5 and 6 we 




R (kpc/h) R (kpc/h) 

Fig. 2. — Density profiles of the most massive ha- 
los at z = 10.5 (full curves). Dashed curves show 
NFW fits. Long-dashed lines in the plot of Halo6 
show power- laws with slopes a = l(NFW) and 
a — 0.2(Ricotti). In spite of large fluctuations at 
large (virial) radii, none of the halos show slopes 
suggested by Ricotti (2002). 

present the NFW parameters Vg and po proper 
units. The name of the halos are given in the last 
column. 
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We impose the following "equilibrium" criteria 
to select a halo for further analysis: (1) The halo 
should not be within the virial radius of a more 
massive halo; i.e. the halo is not a subhalo. (2) 
The halo should not have a subhalo with a mass 
larger than one tenth the mass of the halo. Finally, 
(3) wc also use a goodness of fit criterion to reject 
halos with Df > 0.004, where the accuracy of the 
at Df is defined as follows. 

We compute the density, p{r), in spherical 
shells whose radii increase as the square of the bin 
number, n. Under this scheme, the logarithmic ra- 
dial width, Alogr, is not constant but decreases 
with n, and for high values of n A log r goes as 2/n. 
This binning improves the often used constant log- 
arithmic binning. We use the NFW profile, Pnfw, 
to fit our halo density profiles 



Pnpw(^) — 



Pa 



(1) 



The best fit and thus parameters po and are 
obtained by minimizing x^. 



N 



n=l 



logp(r„) - log 



We define w„ = l/cr^ to be 



Ar 



2n + 7 
(n + 3)2 



(2) 



(3) 



to compensate for the increasing density of bins 
with radius. In equation (3), Ar is the radial 
width for our binning scheme and n is the number 
of the bin corresponding to radius r„. Notice that 
by definition uj is constant for a binning constant 
in log(r). We estimate the accuracy of the fit with 
the parameter D = x^/P^ where F = ^a;„ and 
the sum runs over the number of bins N. This 
parameter reduces to x^/^ (with cr„ = 1) when 
u> is constant and it is a better estimation of the 
goodness of the fit than, for example, a maximum 
deviation because average out possible large fluc- 
tuations in some bins. We set D = Df = 0.004 
as our criterion for the quality of the fit. Halos 
with Df < 0.004 that satisfy at the same time 
our minimum requirements for equilibrium (crite- 
ria 1 and 2) comprise 65% of the sample. In fact, 
in the distribution of D there is a strong drop in 
the number of halos above Df. It seems that Df 
separates relaxed from unrelaxed halos. 




0.1 1 

r (kpc/h) 



10 0.1 



r (kpc/h) 



Fig. 3. — Density and circular velocity profiles of 
the four most massive halos at z = 3.3 that satisfy 
our criteria of "equilibrium" halos. The densities 
of the first (Al, left panel) and third (A3, right 
panel) halos have been multiplied by ten to pre- 
vent overlapping. The corresponding circular ve- 
locities (lower panel) were multiplied by a factor 
of two. The innermost bin is chosen to contain at 
least 200 particles. Vertical lines at the bottom of 
the density panels mark the corresponding virial 
radii (AI, A3 - dotted Unes, A2, A4 soUd lines) 



Figure 3 shows the density and circular ve- 
locity, Vc, profiles for the four most massive ha- 
los of the 1 h~^Mpc box that satisfy the crite- 
ria of "well-behaved" or "equilibrium" halos (Al 
- A4). Dotted and solid lines are the NFW best 
fits. NFW circular velocity profiles are computed 
using {pojfs) parameters obtained from the den- 
sity fits. The NFW profile describes quite well 
both the density and the circular velocity profiles 
of these halos. At redshift 2: 3 the majority - 
about 65% - of dwarf halos with more than 5000 
particles are in "equilibrium" . 

Figure 4 shows the inner structure of the most 
massive halo Al (density, circular velocity and 
3D velocity dispersion profiles) for epochs from 
z = ll.bto z = 2.3. The virial mass of the halo in- 
creases 71 times during the period. Yet, the inner 
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density profile (radii less than w 100 pc) changes 
very little. Thus, the cusp of the halo was already 
in place at very high redshift and since then shows 
little evolution. The central 0.5 kpc region shows 
little evolution since redshift z = 5. This does not 
mean that particles, which at redshift 2; « 10 are 
in the core stay there at later times. Comparison 
of the rms velocities (middle panel) with the cir- 
cular velocities (top panel) clearly indicates that 
particles in the core move with random velocities 
that significantly exceed the circular velocity. This 
implies that the core itself is not self-bound and 
a large fraction of particles simply passes through 
the cusp. This is very important for understand- 
ing the evolution of halo concentration. 




0.0001 



0.05 0.1 



0.5 1 5 10 

Radius (kpc/h) 



Fig. 4. — Circular velocity, 3D velocity dispersion, 
and density profiles of the most massive halo Al 
at different redshifs. Dot dashed curves on the 
top two panels are for z = 3.3. For clarity we 
do not show the density profile for this redshift. 
Vertical lines in the bottom panel mark the virial 
radii. Most of the changes happen at outer radii of 
the halo while there is very little evolution in the 
center. The thick full curve in the bottom panel 
shows the NFW fit for z = 2.3. 

For analysis of halo concentrations we make fits 
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Fig. 5. — The dependence of concentration {cmr, 
top panel), the scale radius (r^, middle panel) and 
Ps = Po/4 {ps, bottom panel) on the halo virial 
mass. Solid squares are halos in the 1 /i~^Mpc 
box at z = 3.3. The rest of the points are for ha- 
los in the 25 h~^Mpc simulation (open triangles 
and four-point stars) and the 80 h~^Mpc simula- 
tion (open squares) at z = 3. The solid and dashed 
curves represent the predictions of the models by 
Bullock et al. (2001a) and Eke, Navarro, & Stein- 
metz (2001), respectively, while the dotted lines 
are the la scatter of halos from the 1 /i~^Mpc box, 
A{logCvir) = 0.11. 



for halos with more than 5000 particles in all our 
simulations. At 2: = 3.3 in the 1 /i~^Mpc box 
simulation there are 72 of these halos of which 
40 satisfy the equilibrium criteria. We also in- 
clude seven more halos that do not satisfy the 
second criterion (have large companion), but are 
still well fitted by the NFW profile. Halos se- 
lected in the other two simulations also have more 
than 5000 particles each and the deviations of the 
measured density profile from the NFW fit are 
less than 30%; \p{ri) - /9NFw(n)| / Pt^Fwin) < 0.3, 
where goes from the first radial bin to r^ir- Ha- 
los that are too disturbed hardly satisfy this latter 
"equilibrium" requirement (Jing 2000). Figure 5 
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shows concentrations Cyir = rmr/fs, proper scale 
radii, r^, and proper densities at the scale radii 
Pa = Po/4: for halos of different virial masses Mvir- 
The solid and dashed curves arc the predictions of 
the concentration at z = 3.0 based on the mod- 
els of Bullock et al. (2001a) and Eke, Navarro, 
& Stcinmetz (2001), respectively, for the appro- 
priate cosmological model with erg = 0.9. The 
points from the 1 h~^Mpc box are slightly scaled 
up to take into account the different normalization 
{as = 0.75) and the slightly different identifica- 
tion redshift {z = 3.3). For this latter, we use the 
1 + z dependence of c-uir while for the former we 
multiply Cyir by 0.9/0.75, which is a good approx- 
imation because dwarf halos form very early when 
the growth factor is proportional to a.^ We use 
the halos from the 1 /i~^Mpc box to compute the 
standard deviation of the log(c) distribution and 
find A{\ogCyir) = 0.11, which is a factor 1.6 lower 
than that quoted by Bullock et al. (2001a). If 
we use the whole sample of halos only under the 
restriction that they must have more than 5000 
particles, the scatter increases to 0.14. So, at least 
part if not all (see §5), of our smaller scatter is due 
to the fact that we are using relaxed halos. 

In summary, dwarf halos seem to be in many 
respects a down-scaled copy of more massive ha- 
los. For example, they all have cusps regardless of 
their equilibrium status. Interestingly, they have 
a median concentration that agree with the model 
of Bullock et al. (2001a). We extended the Cyir 
versus Myir diagram to high masses and showed 
how the Bullock et al. model continue to predict 
the correct concentration. 



2003). This function can be approximated by a 
power-law with power index —2.75 (Klypin et al. 
1999b) (see dot-dashed line in Figure 6). 
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Fig. 6. — The cumulative number of subhalos with 
given circular velocity for halos Al (open circles), 
Bl (open triangles), and CI (solid circles). At 
z ~ 3 each halo has more than 2x10^ particles and 
is resolved down to radius 0.01rvir(the first radial 
bin with > 200 particles). The subhalo maximum 
circular velocities are in units of the virial velocity 
of the parent halo. The dot-dashed line shows the 
power law with the slope -2.75 found in previous 
simulations. 



3.2. Substructure mass function 

The cumulative velocity function of subhalos 

within halos Al (open circles), Bl (open trian- 
gles), and CI (solid circles) are compared in Figure 
6. The peak velocities of subhalos are measured in 
units of the virial velocity, defined as the circular 
velocity at the virial radius, of their parent halo. 
The mass of the halos covered almost four orders of 
magnitude yet the subhalo velocity distributions 
are similar (Moore et al. 2001; De Lucia et al. 

^Bullock et al. (2001a) define the epoch of collapse as the 
epoch at which the non-lineaj: mass, M«, equals a fixed 
fra<;tion F of the halo mass, Mvir- For the ACDM model, 
Bullock et al. suggest to use F = 0.01. 



4. Statistics of halos 
4.1. Halo mass function 

The mass function of halos N{M) is defined so 
that N{M)dMdV represents the number of ha- 
los with masses between M and M + dM in the 
volume dV. This is a differential mass function. 
In order to find isolated halos we use the friends- 
of-friends group finding algorithm (Davis et al. 
1985) with a linking length of 0.2 (or in physi- 
cal units 0.78 h~^kpc). Halos are identified as the 
set of particles mutually linked. The mass of the 
halo is simply the sum of its constituent dark mat- 
ter particles. Tests (e.g., Governato et al. 1999) 
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indicate that at least 30 particles arc needed to 
take into account properly the effect of mass res- 
olution. In this paper we analyze only halos with 
more than 40 particles. 

We construct the dwarf halo mass function by 
grouping the halos according to their mass in bins 
of Alog(M) = 0.1. In order to avoid too large 
statistical fluctuations, we do not use the few most 
massive halos. Error bars for data in individual 
bins are estimated using Poisson statistics. 




M( Mg h-') 

Fig. 7. — The halo mass function from the 1 
/i~^Mpc box (open squares). Error bars are com- 
puted assuming the Poisson statistics. The solid 
line is the analytical prediction from the Sheth- 
Tormen formalism. 

Figure 7 shows the mass function (open squares) 
of our dwarf halos from the 1 h~^Mpc box. The 
solid line is the analytical prediction of the mass 
function from the Sheth-Tormen formalism. No- 
tice the excellent agreement between the predic- 
tion of the Sheth-Tormen model and our numerical 
results. This model has been tested under diverse 
conditions: redshifts, cosmological parameters, 
halo masses (> lO^°/i~^M0), etc., and have passed 
all tests (but see Reed et al. 2003). Here we see 
that it successfully predicts the mass function of 
dwarf halos, albeit at high redshift. In summary, 
the Sheth-Tormen formalism can be safely applied 
in different cosmological models and for about 
ten orders of magnitude in mass. With the cor- 
responding correction factors the Sheth-Tormen 
formalism can be successfully applied also to pre- 



dict the number of dwarfs in voids (Gottlober et 
al. 2003). 

4.2. Spin parameter distribution 

To study the distribution of the spin parameter, 
p(A), we use halos which do not reside inside a 
larger halo and have more than 500 particles. The 
spin parameter A is defined by 



A; 



J\E 



1/2 



GM5/2 ' 



(4) 



where J is the magnitude of the angular momen- 
tum, E is the total energy, and M is the mass of 
the halo. We assume that halos are in virial equi- 
librium and thus E = —K, where K is the kinetic 
energy inside rvir- N-body simulations have shown 
that p(A) is well described by a lognormal distri- 
bution (see Vitvitska et al. (2002) and references 
therein) 



p{X)d\ = 



1 



exp 



In' (A/Ao) 



2(jx 



Figure 8 compares the probability distributions of 
the spin parameter in simulations A and D (Ta- 
ble 1) and shows the lognormal fits. 




Fig. 8. — Spin parameter distribution (histogram) 
for halos with more than 500 particles drawn from 
a 60 h~^Mpc box simulation at z = (left panel) 
and from our 1 /i~^Mpc box at z = 3.3 (right 
panel) . Curves on each panel are lognormal best- 
fits to the data. The parameters of the fit {a\, Aq) 
are also shown in panels. 

Figure 8 shows that the distribution of the spin 
parameter of dwarf halos also follows a lognormal 
law. In fact, it turns out that dwarf halos in simu- 
lations A are better described by equation (5) than 
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the more massive halos in the simulation D: their 
distribution has a lower for the same number 
of bins. According to the values of the parameters 
of the fit, it would appear that the distributions 
are different, but within the errors they are simi- 
lar. What happens is that p{X) is very sensitive to 
(cta, Aq). In general, our results confirm the known 
independence of the spin parameter on mass (Lem- 
son & Kauffmann 1999) and epoch (e.g., Vitvitska 
et al. 2002). 

5. Discussion and conclusions 

We find that the density profiles of relaxed 
dwarf halos can be well fitted by the NFW pro- 
file, which has an inner density slope a = —1. 
This result is in conflict with results of Ricotti 
(2002), who finds a « —0.5 for dwarf halos in his 
1 h~^Mpc box simulation. In a attempt to find the 
cause of our discrepancy, we studied density pro- 
files at the same redshift as Ricotti did {z — 10.5) 
and reproduced Ricotti's fitting procedure. The 
conclusion is that even without any fits none of 
the profiles are close to the shallow profiles sug- 
gested by Ricotti. On the other hand, Navarro et 
al. (2003), in a recent study of mass profiles of 
ACDM halos at z = that span five dec;adc;s in 
halo mass, show that inner density slopes of dwarf 
halos are steeper than —1 all the way down to the 
resolution limit. From their Tables 2 and 3 we 
find dwarf halo concentrations from Cyir = 16.4 
for halo Dl to Cmr = 22.0 for halo D4. These con- 
centrations are below to our extrapolated value of 
the median concentration we find for our dwarf ha- 
los, but once the mass dependence of the concen- 
tration is considered they turn out to be similar. 
Our results thus agree with the study of Navarro 
et al. on dwarf halos not only in the inner shape 
of density profiles but also on concentrations. 

It is impressive to see how well the model of 
Bullock et al. (2001a) reproduces over five or- 
ders of magnitude the measured concentration as 
a function of halo mass. At the same time. Bul- 
lock et al. overpredict the deviations from the 
mean value. We compute the median and stan- 
dard deviation of the c-uir distribution of our sam- 
ple of relaxed dwarf halos and find a spread which 
disagrees with that reported in Bullock et al. by 
about a factor 1.6 in Alogc^,;^- It appears that the 
printed version of Bullock et al. has a typo and the 



actual value for the spread of concentration in Bul- 
lock et al. should have been A log c„ir = 0.14 (Bul- 
lock, private communication). This is still larger 
than what we found: AlogCyir = 0.11. Our la 
scatter is close to Alogc^j^ = 0.10 measured by 
Klypin et al. (2003). We found that, except at the 
high-mass end, Eke's et al. model lie below Bul- 
lock's et al. For instance, at Myir = lO^/i~^M0 
Bullock's et al. model predicts Cyir = 8.8 as com- 
pared with 5.5 of Eke's et al. model. 

If we extrapolate the concentration of our dwarf 
halos to the present using the 1+z growth law of 
the Bullock et al. model, the median would be 
about 30. This is higher; for example, than the 
value C200 ^ 14 (or Cvir ~ 18, no mass loss case) 
used by Hayashi et al. (2003). A shift up in C200 
in Hayashi et al. fitting procedure would reduce 
the lower limit in Wmax; that is, a lower t^max would 
be required to match the observed stellar velocity 
dispersion of Carina and Draco. In any case, it 
is not clear to what extent Hayashi et al. results 
can be applied to the observed satellites when sub- 
structure is constrained to evolve in a static poten- 
tial. Stoehr et al. (2002), on the other hand, are 
able to explain the observed central velocity dis- 
persions of the dSph satellites of the Milky Way 
for the 20 most massive subhalos of a Milky Way- 
sized galaxy halo. Subhalos are resolved with hun- 
dreds or at most with a few thousands of particles. 
These still small numbers of particles are proba- 
bly not enough to draw yet firm conclusions. In 
the light of our results, it seems that subhalos re- 
solved with many more particles than subhalos in 
present-day cosmological simulations will be more 
concentrated. However, it is not clear how this ef- 
fect will affect the conclusion that observed dSph 
galaxies inhabit the largest subhalos. 

In our analysis of the structure of dwarf ha- 
los, we also compute the abundance of subhalos. 
We measure the amount of substructure in halos 
Al (dwarf halo), Bl (Magellanic Cloud-sized halo) 
and CI (group-size halo) through the normalized 
subhalo Umax function (umax is measured in units 
of the virial velocity of the par(uit halo) and find 
a mass-scale independent result: the curves are 
similar and well approximated by the power-law 
n(> V) oc (Klypin et al. 1999b). This 

is in agreement with the paper by Moore et al. 
(2001) (see also De Lucia et al. (2003)). 

We compute the dwarf halo mass function and 
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find an excellent agreement with the model of 
Sheth & Tormen (1999). The model has been 
tested at different redshifts, for a variety of cos- 
mologies, and for halo masses above ^ 10^'^h~^MQ 
(Sheth & Tormen 1999; Jenkins et al. 2001; 
Reed et al. 2003). It provides an excellent match 
to the data. In this paper, we tested it at much 
lower masses, 10^ — 10^h~^MQ, and found also an 
excellent agreement with the data. It is pleasant 
to know that this approximation is able to match 
data for ten orders of magnitude in mass. 

The spin parameter distribution p(A) is mea- 
sured for our sample of dwarf halos (halos with 
a virial mass between 10^ and 10^ h~^MQ). It 
is well fitted by a lognormal distribution with 
Ao = 0.042 and a\ = 0.63. These parameters 
roughly agree with those estimated by Vitvitska 
et al. (2002). Our results corroborate the known 
independence of the spin parameter on mass (Lem- 
son & Kauffmann 1999) and epoch (e.g., Vitvitska 
et al. 2002). 
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